Dynamics of axial separation in long rotating drums 
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We propose a continuum description for the axial separation of granular materials in a long rotating 
drum. The model, operating with two local variables, concentration difference and the dynamic angle 
of repose, describes both initial transient traveling wave dynamics and long-term segregation of the 
binary mixture. Segregation proceeds through ultra-slow logarithmic coarsening. 

PACS: 46.10,+z, 47.54.+r, 47.35.+i 



The collective dynamics of granular materials recently 
have been attracting much interest (l]-|(| . Intrinsic dissi- 
pative nature of interaction among macroscopic particles 
sets granular matter apart from familiar gaseous, liquid, 
or solid states. One of the most fascinating features of 
heterogeneous granular materials is their ability to segre- 
gate under external agitation instead of mixing as would 
be expected from thermodynamics. Essentially any vari- 
ation in mechanical properties of particles (size, shape, 
density, surface roughness) may lead to their segrega- 
tion. Long rotating drums partially filled with a mixture 
of grain sizes exhibit both radial and axial size segrega- 
tion In radial segregation grains of one type (for 
grains of different sizes, the smaller ones) build up a core 
near the axis of rotation. Axial segregation often follows 
radial segregation, with the mixture of grains separating 
into bands according to size arranged along the axis of the 
drum. Axial segregation leads to either a stable array of 
concentration bands, or to the complete segregation j7j|. 

Granular flow in a rotating drum is different from con- 
ventional fluid motion. In the bulk, particles perform 
solid body rotation around the drum axis until they reach 
the free surface. Then they slide down within a thin near- 
surface layer ||. The surface has a 5-curved shape and 
its average slope is determined by the dynamic angle of 
repose. Since there is almost no inter-particle motion in 
the bulk, the segregation predominantly occurs within 
the fluidized near-surface layer. The radial segregation 
occurs during first few revolution of the drum. For long 
drums (length much larger than radius), along with ra- 
dial segregation, axial segregation occurs at much later 
stages (after hundreds of revolutions). Recent experi- 
ments IIHD revea l e d interesting features of axial segre- 
gation. At early stages, small-scale perturbations travel 
across the drum in both directions (it was obvious from 
the dynamics of pre-segregated mixtures (6|), while at 
later times more long-scale static perturbations take over 
and and lead to quasi-stationary bands of separated mix- 
ture. Bands of segregated materials interact at a very 
long time scale and exhibit logarithmic coarsening ||. 



This process can be accelerated in a drum of a helicoidal 
shape. The bands can be locked in a drum with the ra- 
dius modulated along the axis |J . 

Most of the theoretical models of segregation agree in 
that the underlying reason for segregation is the sensi- 
tive dependence of the surface slope or shape on the rel- 
ative concentration of different particles in the mixture 
||. In Ref. H a simple theory was proposed of segrega- 
tion in thin surface flow driven by the local slope, com- 
bined with conservation of particles. For monodisperse 
material, the model recovers the S'-shape of the free sur- 
face. For a binary mixture the model yields a nonlinear 
diffusion equation for the relative concentration of the 
ingredients along horizontal axis. Axial segregation oc- 
curs when the diffusion coefficient turns negative. This 
model yields a significant insight into the nature of the 
instability leading to the segregation, however since it is 
based on a first-order diffusion equation, it fails to de- 
scribe the traveling waves observed at the early stages of 
axial segregation B. 

In this Letter we propose a continuum model which 
describes consistently the early phase of segregation with 
traveling waves as well as the later stage of segregation 
characterized by slow merging of bands of different parti- 
cles. Our model predicts slow (logarithmic) coarsening of 
the segregated state. The dynamics of segregation shows 
striking similarity with the experiments of Ref. j(| . 

Let us consider a mixture of two sorts of particles, A 
and B, of which A(B) corresponds to particles with larger 
(smaller) static repose angle. Our model operates with 
two continuum variables: the relative concentration of 
particles c = (ca — c b)/(c), and local dynamic repose 
angle 0. Here ca,b are local partial concentrations of 
particles, and (c) = {ca + cb) is an average over whole 
system total concentration. We assume that c and are 
functions of longitudinal coordinate x and time t. In 
fact, concentration also depends on the radial coordinate, 
however in this study we are concerned with axial segre- 
gation, and will operate with quantities integrated over 
cross-sections of the long drum. 
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The first equation represents conservation of the rela- 
tive concentration c in the binary mixture: 

d t c=-d x (-Dd x c + g(c)d x 9) (1) 

The first term describes diffusion flux (mixing), and the 
second term describes differential flux of particles due to 
the gradient of dynamic repose angle. Similar term with 
specific function g(c) = Gq(1 — c 2 ) (in our notation) was 
first derived in Ref. || from the condition of detailed flux 
balance. The proportionality constant Go depends on the 
physical properties of the grains (see ||). For simplicity, 
the constant Go can be eliminated by the scaling of dis- 
tance x — > x/\/Gq. The sign + before this term means 
that the particles with the larger static repose angle are 
driven towards greater dynamic repose angle. As will be 
shown below, this differential flux gives rise to the seg- 
regation instability. Since this segregation flux vanishes 
with g(c) \c\ — > 1 (which correspond to pure A or B 
states), it provides a natural saturation mechanism for 
the segregation instability. 

The second equation describes the dynamics of 9 

d t e = a (n-e + /(c)) + D B d xx e + jd xx c. (2) 

Here fl is the angular velocity of the drum rotation, and 
/(c) is the static angle of repose which depends on the 
relative concentration. According to our definition of 
c, /(c) is an increasing function of relative concentra- 
tion Q . Since the angle of repose as a function of the 
concentration c changes typically in a small range, we 
can approximate the function /(c) by linear dependence 
/(c) = F + foe. The constant F can be eliminated by 
the substitution 9 — ► 9 — F. First term in the r.h.s. of 
Eq.(|]) describes the local dynamics of the repose angle 
(D, increases the angle, and —9 + /(c) describes the equi- 
librating effect of the surface flow), and the term Dgd xx 9 
describes axial diffusive relaxation. The last term, jd X xC, 
represents the lowest-order non-local contribution from 
the inhomogeneous distribution of c (the first deriva- 
tive d x c cannot be present due to reflection symmetry 
x — > — x). As we will see later, this term gives rise to the 
transient oscillatory dynamics of the binary mixture. 

The stationary uniform state of the system is c = 
c o; #o = fi + /oco where Co is determined by initial condi- 
tions. 9q has a meaning of the stationary dynamic repose 
angle, which in the limit O — > approaches the static re- 
pose angle 0. 

Let us consider the stability of the uniform state, 

c = co + Ce xt+lkx , 

9 = 9 + $e xt+lkx (3) 

To linearize the system, we need to expand the function 
g(c) near Co. The stability properties depend on the val- 
ues of /o and go = g(co). The eigenvalues are found 
from the following equation, 



(A + Dk 2 ){\ + Dgk 2 + a) - g (a * /Ofc 2 - 7 fc 4 ) = 0. (4) 

At fc — > one finds Ai = —a — (Dg + gofo)k 2 ; A2 = 
(50/0 — D)k 2 . Asymptotic expansion at k — » 00 yields 

Ai, 2 = \ (-Dg -D± ^{Dg - Df - 4507) k 2 . (5) 

It is easy to see that if gofo > cxD, long-wave perturba- 
tions are unstable (A 2 > 0), and if go7 > (Dg — D) 2 /4, 
short-wave perturbations oscillate and decay (eigenval- 
ues X12 are complex conjugate with negative real part). 
The full dispersion curve A(fc 2 ) for particular values of 
parameters is plotted in Fig. [I]. This curve is consis- 
tent with the measurements in Ref. E| . Indeed, at small 
wavelengths the perturbations travel, and their phase ve- 
locity v p h — ImA/fc grows linearly with k. At small fc, 
frequency and phase velocity turn into zero and remain 
zero for all long- wave perturbations k < fc*. Perturba- 
tions in the range < fc < fc c where k c = (ce(go/o — 
D)/(gol + DDg)) 1 / 2 < fc*, grow exponentially and there 
is an optimal wavenumber of the fastest growing pertur- 
bations fco- Oscillating perturbations with fc > fc* always 
decay, which seems to be in agreement with experimen- 
tal data, although direct experimental measurement of 
ReA(fc) is lacking. 

We performed numerical simulations of Eqs.(0j^) using 
pseudo-spectral split-step method with periodic bound- 
ary conditions. We used up to 512 mesh points in our 
numerical procedure. The following form of nonlinear 
functions in Eqs. (HJ|) was implemented: g[c) = 1 — c 2 
and /(c) = f c |10| . The dynamics of the initially pre- 
separated state with wavenumber fc > fc* in a system with 
size L = 60 is shown in Fig. |2|. In agreement with experi- 
ments II , short-wave initial perturbation produce decay- 
ing standing waves, which later are replaced by quasi- 
stationary bands (Fig. ||a). The bands are separated 
by sharp interfaces which are very weakly attracted to 
each other. In fact, in simulations with parameters cor- 
responding to Fig.|2|a we were not able to detect interface 
merging at all in a reasonable simulation time. However 
at higher rates of diffusion and dissipation, corresponding 
to different experimental conditions (rotation frequency, 
filling factor, etc), the interaction becomes more signif- 
icant, and it leads to band merging and overall pattern 
coarsening (see Fig. ||d). In Fig. [| we present a number 
of bands as a function of time for this run. 

Fronts separating bands of different grains, can be 
found as stationary solutions of Eqs. (0) , (§) . In an infinite 
system one finds from stationary Eq.([I|) 9 = 9q+ DG(c), 
where G(c) = /[g(c)] _1 dc = — ~ log and 9 is an inte- 
gration constant. Plugging this expression in Eq.(||), we 
obtain the second order differential equation for c (for a 
symmetric solution one chooses 9q = Q), 

d \( DDg \ del „ aD , 1-c 
Tx[h + —)dx\ +af ° C+ - l °Z—= - (6) 
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It is easy to see that this equation possesses an interface 
solution. The asymptotic behavior of this solution can 
be found in the limit D <C /o, when the states on both 
sides of the interface are well segregated (\c(x — > oo)| — > 
±1). In this limit, far away from the interface 1 — |c| oc 

cxp(-x/d ) where d a = yj ' + ^ exp(-^). As seen 
from this formula, the characteristic front width vanishes 
as D, Dg — > 0. This result could be anticipated, as in the 
absence of diffusion nonlincarity g(c) drives the system 
towards complete segregation. 

Slow logarithmic coarsening of the segregated state can 
be understood in terms of the weak interaction of fronts. 
Since the asymptotic field of the front approaches equi- 
librium value of the concentration exponentially fast, we 
expect in general exponentially weak interaction between 
the neighboring fronts leading to logarithmic times for 
the front annihilation. 

The problem of front interaction has been recently con- 
sidered |l2j for the Cahn-Hilliard model which arises in 
relaxation phase ordering kinetics pif . As in the Cahn- 
Hilliard model, the order parameter (here concentration 
c) is a conserved quantity, therefore front interaction 
must conform to a global constraint. In fact, in a partic- 
ular case g(c) = 1, D = 0, Dg = our model (0),(§) can 
be reduced to a single equation for the relative concen- 
tration 

d tt c + d t c = -d xx (f(c) + d xx c) (7) 

which differs from the Cahn-Hilliard model only by the 
first term in the l.h.s. This term describes non-potential 
effects, including traveling waves discussed above. For 
the slow process of the band relaxation in the long-time 
limit this term can be dropped. Front solutions within 
the Cahn-Hilliard equation exist for iV-shaped functions 
/(c), for example, a cubic polynomial /(c) = 2(c— c 3 ). 
In this case stationary interface solution is found analyt- 
ically as c(x) = Co + tanh(x — xq), and interaction of the 
interfaces can be analyzed. This analysis [j^], predicts 
that a single band of positive c in the sea of negative c 
(two fronts) can either annihilate or reach a stationary 
width depending on the initial distance between fronts. 
If a number of fronts is greater than 2, front interaction 
leads to their annihilation and pattern coarsening. The 
interaction between fronts due to their tail overlapping 
is exponentially weak. The rate of front interaction can 
be determined using multiple scale analysis. In the sim- 
plest case of three equidistant fronts, two outer fronts 
collide with the central front and disappear after time 
T oc dexp(d/do) where d is the initial distance between 
fronts. Thus, for multi-band structure, the number of 
fronts N (proportional to the inverse average distance 
between fronts) decreases approximately as a logarithmic 
function of time (N = 1/d ~ 1/ (const + do logT). Since 
interfaces in our more complicated model have similar 
exponentially small tails, we anticipate a similar loga- 



rithmic law of band coarsening. This dependence indeed 
agrees with our numerical simulations (see Fig. |J). 

In conclusion, we proposed a simple continuum model 
for axial segregation of binary granular mixtures in long 
rotating drums. The model operates with two local dy- 
namical variables, relative concentration of two compo- 
nents and dynamic repose angle. The dynamics of our 
model shows qualitative similarity with the experimen- 
tal observations of initial transients and long-term seg- 
regation dynamics |^]J^^|- It captures both initial tran- 
sient traveling waves and subsequent onset band struc- 
ture. The dispersion relation for slightly perturbed uni- 
form state (f|) qualitatively agrees well with observations 
|^|, and can serve for fitting the model parameters. The 
logarithmic coarsening of the quasi-static band structure 
which follows from our model, deserves experimental ver- 
ification. This coarsening is typical for systems with ex- 
ponentially weak attractive interaction among defects or 
interfaces, as in the phase ordering kinetics described by 
the Cahn-Hilliard model. Our simulations also showed 
that the model qualitatively reproduces more compli- 
cated phenomenology of the separation process reported 
in Ref. In particular, periodic modulation of the 

drum radius, modeled in our approach by periodic varia- 
tion of f2, leads to band locking. Breaking of the x — > —x 
symmetry by term d x 9, introduced in r.h.s. of Eq. (Q), 
results in complete segregation, similarly to the dynamics 
of grains in the drum with hclicoidal shape. 
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FIG. 1. Dispersion relation X(k) for the small per- 
turbation of the uniform state cq = at /o = 40, 
D = 0.05,A> = 0.1, a = 0.025, 7 = 1. Perturbations are 
unstable at k < k c — 0.994 and oscillate (and decay) at 
k > fc* = 1. 




FIG. 2. Space-time diagrams of the evolution of the 
pre-separated state (c(x,0) = Co + Ci cos(kox)) with 
co = 0., Ci — 0.95, and fco = 1.79 > k*);a ini- 
tial transient, At times t < 15 the initial perturbations 
excite decaying standing wave (superposition of left- and 
right-traveling waves), and at larger times t > 15, ape- 
riodic segregated bands emerge. Parameters of the model 
are: D = 0.05,7 = 1,Q = 0.025, / = 40, L = 60 and 
De — 0.1. b, Space-time diagram for long-time evolution, 
band merging and coarsening during long-time evolution 
(0 < t < 10000) at higher diffusion constants, parameters 
D = 0.8, D e = 0.5, 7 = 1, a = 0.5, fo=2, L = 140. 
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FIG. 3. Number of fronts N as a function of time (di- 
amonds) and its fit by a function iV = 70/(ln£ — 2.5) 
(long-dashed line). Parameters correspond to Fig. 2b 
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